A phase field approach for direct calculation of interfacial free energy of solid-gas and solid-liquid interfaces in Lennard-Jones systems

Direct calculation of solid-liquid interfacial free energy for molecular systems: TIP4P ice-water interface.

By extending the cleaving method to molecular systems, we perform direct calculations of the ice Ih-water interfacial free energy for the TIP4P model. The values for the basal, prism, and {112[over]0} faces are 23.3+/-0.8 mJ m{-2}, 23.6+/-1.0 mJ m{-2}, and 24.7+/-0.8 mJ m{-2}, respectively. The closeness of these values implies a minimal role of thermodynamic factors in the anisotropic growth o...

Relaxation of surface tension in the liquid-solid interfaces of Lennard-Jones liquids.

We have established the surface tension relaxation time in the liquid-solid interfaces of Lennard-Jones (LJ) liquids by means of direct measurements in molecular dynamics (MD) simulations. The main result is that the relaxation time is found to be almost independent of the molecular structures and viscosity of the liquids (at 70-fold change) used in our study and lies in such a range that in sl...

The solid-liquid interfacial free energy out of equilibrium

The properties of the interface between solid and melt are key to solidification and melting, as the interfacial free energy introduces a kinetic barrier to phase transitions. This makes solidification happen below the melting temperature, in out-of-equilibrium conditions at which the interfacial free energy is ill-defined. Here we draw a connection between the atomistic description of a diffus...

Solid-liquid interfacial free energy out of equilibrium

Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation.

Computing phase diagrams of model systems is an essential part of computational condensed matter physics. In this paper, we discuss in detail the interface pinning (IP) method for calculation of the Gibbs free energy difference between a solid and a liquid. This is done in a single equilibrium simulation by applying a harmonic field that biases the system towards two-phase configurations. The G...